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Abstract 


The  electrical  and  rheological  properties  of  two  state-of-the-art  electro-rheological 
(ER)  fluids  have  been  evaluated  from  0  to  60  °C,  using  a  unique  constant-stress 
viscometer.  The  measurements  were  performed  at  the  facilities  of  ER  Fluid 
Developments,  Ltd.  (ERFD),  a  small  British  company  located  in  Sheffield,  United 
Kingdom,  that  specializes  in  developing  and  marketing  ER  fluid-based  technology.  In 
addition  to  these  measurements,  the  response  time  of  the  fluids  was  determined  using 
another  instrument  located  at  the  nearby  University  of  Sheffield. 

The  constant-stress  viscometer,  or  dynamic  rig,  developed  and  owned  by  ERFD, 
permits  the  yield  stress,  plastic  viscosity,  and  direct  current  (DC)  current  density  of  an 
ER  fluid  to  be  determined  as  a  function  of  electric  field  and  temperature.  The  yield 
stress,  plastic  viscosity,  current  density,  and  response  time  of  each  fluid  were  found  to 
strongly  depend  on  temperature.  The  first  fluid  appears  to  be  stronger  (i.e.,  has  a  higher 
yield  stress),  with  lower  current  densities  at  lower  temperatures  up  to  40  °C.  In  that 
temperature  range,  the  second  fluid  is  somewhat  weaker,  susceptible  to  electrical 
breakdown  at  or  above  2.0  kVmm  1,  and  generally  draws  more  current.  At  higher 
temperatures,  however,  the  current  density  of  the  first  fluid  dramatically  increases  to 
exceed  the  desired  level  of  10  jiAcm'2.  The  corresponding  current  density  increase  of 
the  second  fluid  is  significantly  smaller.  In  light  of  the  differences  between  the 
measured  properties,  the  results  and  their  implications  associated  with  the  intended  use 
of  the  fluids  are  discussed.  A  general  description  of  ER  fluid  behavior,  as  well  as  the 
theories  of  operation  of  the  instruments,  are  also  described. 
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1.  Introduction 


Several  Army  laboratories,  with  the  lead  of  the  Fire  Support  Armaments  Center, 
U.S.  Army  Armaments  Research,  Development,  and  Engineering  Center 
(FSAC/ARDEC),  Rock  Island  Arsenal  (RIA),  U.S.  Army  Research  Laboratory 
(ARL),  and  General  Atomics  Corporation  (GA)  are  currently  developing  a 
large-caliber,  155-mm  Light  Artillery  Future  Direct  Support  Weapon  System 
(FDSWS)  for  Army  After  Next  forces.  Since  its  inception,  the  FDSWS  has  been 
renamed  to  the  Advanced  Technology  Lightweight  Artillery  System  or  ATLAS. 
Relying  on  advanced  materials  and  novel  recoil  abatement  technologies,  ATLAS 
is  intended  to  provide  lethality  and  firepower  without  compromising  accuracy, 
weight,  or  mobility. 

One  of  the  key  elements  of  ATLAS  is  its  use  of  electro-rheological  (ER)  fluids  to 
control  the  recoil  action  of  the  gun  barrel,  wherein  the  impulse  energy  imparted 
by  the  recoil  stroke  is  managed  and  dissipated  over  an  expanded  time  span.  A 
benefit  to  using  ER  fluids  either  directly  or  indirectly  (used  as  the  control  fluid 
for  the  otherwise  nominal  howitzer  hydraulics)  is  that  when  the  peak  stresses  are 
significantly  reduced,  the  need  for  the  massive  undercarriage  and  supports  are 
lessened. 

For  ATLAS,  an  ER  fluid  with  high  activity  (high  shear  strength),  low  no-field 
viscosity,  low  current  density,  and  fast  response  time  over  a  wide  temperature 
range  is  desired.  Additionally,  the  fluid  should  be  nontoxic,  noncorrosive, 
nonsettling,  and  compatible  with  a  variety  of  common  engineering  materials. 
Table  1  summarizes  these  "benchmark,"  minimum  property  values. 


Table  1.  Desired  benchmark  properties  of  candidate  ER  fluids. 


Property 

Target  Value 

Conditions 

Activity  or  Yield  Strength;  xy 

4  kPa 

4  kV/mm;  20  °C 

Direct  Current  (DC)  Density;  j 

<  10  pA/  cm2 

0-40  °C 

Zero-Field  Viscosity;  rjo 

<  120mPa»s 

20  °C 

Response  Time;  tr 

2ms 

0-40 °C 

Nontoxic,  noncorrosive,  and  nonsettling 

Compatible  with  Aluminum,  Titanium,  Stainless  Steel,  Delrin,  etc. 

In  the  past,  ER  fluids  have  been  applied  to  automotive  engine  mounts,  clutches, 
valves,  vibration  and  noise  reduction  control  mechanisms,  and  liquid  crystal 
displays.  However,  there  are  currently  still  no  significant  military  applications  in 
use  [1].  Although  several  fluids  are  now  commercially  available,  they  tend  to  be 
highly  design  specific.  That  is,  due  to  its  highly  specific  properties,  one  fluid 
used  in  one  application  may  not  be  suitable  in  another.  Furthermore,  data 
representing  the  fluid  characteristics  is  often  incomplete,  so  an  accurate  property 
assessment  cannot  be  made.  Occasionally,  such  data  may  not  even  be  available 
from  the  supplier  because  of  a  lack  of  proper  instrumentation  for  its 
determination.  Consequently,  ARL  has  been  tasked  to  develop  an  in-house 
evaluation  capability  to  assess  ER  fluids  on  a  comparative  basis,  as  well  as  create 
a  comprehensive  engineering  property  database  for  ER  fluids. 

Since  the  FY95-FY96  timeframe,  GA  worked  under  contract  with  FSAC  towards 
designing,  constructing,  and  testing  a  lightweight,  full-scale,  ER  fluid-based, 
soft-recoil  arrest  simulator  and  its  associated  control  system  [2,  3].  With  the  soft 
recoil  system  fitted  onto  a  test  bed,  live-fire  testing  was  to  commence  to 
demonstrate  the  feasibility  of  closed-loop,  real-time  control  of  a  howitzer  using 
ER  fluids.  Previously,  GA  and  its  subcontractor,  ER  Fluid  Developments,  Ltd. 
(ERFD),  in  Sheffield,  United  Kingdom,  have  already  built  and  tested  a 
proof-of-principle  subscale  device  for  FSAC  that  has  essentially  demonstrated 
the  feasibility  of  the  ER  fluid  soft-recoil  concept  [3].  ERFD  has  been  developing 
ER  fluids  and  associated  applications  for  over  30  years. 

In  the  fall  of  1996,  a  request-for-quote  (RFQ)  was  placed  by  GA  for  an  ER  fluid 
that  could  meet  the  benchmark  requirements  of  die  ATLAS  program.  Two 
candidate  samples  were  submitted  in  December  1996.  GA  was  to  downselect 
one  of  the  fluids  for  ATLAS  after  a  thorough  property  assessment.  GA  requested 
that  ARL,  acting  as  an  independent,  disinterested  "third"  party,  measure  the  key 
physical  properties  of  the  fluids.  However,  at  that  time,  ARL  still  lacked  the 
proper  instrumentation  and  facilities  for  the  evaluation.  Therefore,  GA 
suggested  and  arranged  for  ARL  to  perform  the  necessary  measurements  by 
working  with  ERFD  in  England. 

Test  requirements,  instructions  for  handling,  and  fluid  samples  were  shipped  to 
ARL  and,  in  February  1997,  the  author  traveled  to  ERFD  to  perform  the 
measurements.  The  fluid  identities  and  their  respective  suppliers  were  known 
only  to  GA.  This  was  crucial  for  providing  a  comparative  but  anonymous  basis 
for  experimentation  while  maintaining  confidentiality  and  impartiality. 

Because  of  time  restrictions,  the  evaluation  effort  primarily  concentrated  on  the 
properties  of  shear  stress,  viscosity,  and  current  density  as  a  function  of  electric 
field  intensity,  shear  rate,  and  temperature.  ERFD's  constant  stress  rheometer,  or 
dynamic  drop  viscometer,  was  used  to  measure  the  current  density,  plastic 
viscosity,  and  shear  strength  of  the  ER  fluids  as  a  function  of  temperature.  In  a 
separate  series  of  experiments,  the  response  time  of  the  fluids  was  measured 
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under  a  variety  of  conditions.  These  latter  tests  were  performed  at  University  of 
Sheffield's  Department  of  Mechanical  and  Process  Engineering,  Smart  Machines, 
Materials,  and  Related  Technologies  (SMMART)  Laboratory.  Similar  to  the 
viscometry  experiments,  the  scope  of  these  tests  was  to  examine  the  effects  of 
shear  rate,  electric  field  intensity,  and  temperature.  This  report  enumerates  the 
results  of  both  of  these  series  of  experiments  and  describes  the  devices  used  for 
their  determination. 


2.  ER  Fluid  Characteristics 


This  section  is  primarily  intended  to  provide  a  basis  for  interpretating  test  results 
by  the  uninitiated.  The  subject  of  ER  fluids  has  been  extensively  covered  over 
the  last  30  years,  and  several  excellent  review  articles  have  appeared  in  the 
literature  [4-9]. 

Currently,  only  a  handful  of  manufacturers  in  the  United  States  and  abroad  have 
the  capability  to  produce  batch-size  quantities  (i.e.,  more  than  a  few  liters  at  a 
time)  of  ER  fluids.  In  the  United  States,  Bayer  Corporation  (Pittsburgh,  PA)  and 
Bridgestone-Firestone  Corporation  (Akron,  OH)  are  still  active,  in  Europe,  DEA 
Mineraloel  AG  (Hamburg,  Germany)  and  ERFD  (Sheffield,  United  Kingdom) 
continue  to  produce  fluids.  Actually,  in  practice,  ERFD  has  made  arrangements 
with  Loctite  Corporation  (Dublin,  Ireland)  for  the  exclusive  license  to 
manufacture  and  sell  larger  quantities  of  ERFD's  fluid  formulations.  In  Japan, 
Asahi  Chemical  Industry  Company  (Fuji-shi,  Shizuoka-ken),  Fujikura,  Ltd. 
(Koto-ku,  Tokyo),  and  Nippon  Shokubai  Corporation  (Tsukuba,  Ibaraki)  lead 
among  other  manufacturers.  Former  fluid  producers  and  their  subsidiaries,  such 
as  Lord  Corporation  (Cary,  NC),  Lubrizol  (Wickliffe,  OH),  RWE  (Wesseling, 
Germany),  and  BASF  (Ludwigshafen,  Germany),  etc.,  have  either  discontinued 
their  research  and  development  efforts,  or  are  continuing  them  on  a  greatly 
scaled-back  level. 

2.1  General  Trends 

ER  fluids  are  fine  particle  suspensions  in  a  base  liquid.  By  applying  an  electric 
field,  they  increase  in  viscosity  to  the  point  of  becoming  almost  solid-like, 
thereby  capable  of  supporting  higher  shear  stress  levels  than  normal  fluids.  The 
composition  of  the  conductive,  polar  particles  can  range  from  flour,  starch, 
crystalline  cellulose,  carbon,  polymers  with  or  without  attached  metal  ions, 
metals,  or  ceramics  (e.g.,  alumina,  silica,  titania,  zeolites,  etc.).  The  concentration 
of  particles  can  vary  from  10  to  50  volume-percent,  while  the  optimum  particle 
size  seems  to  be  between  40  run  and  50  pm.  At  one  extreme,  interference  from 
the  Brownian  motion  of  ultrafine  particles  limits  the  effectiveness,  but  at  the 
other  extreme,  large  particles  suffer  from  sedimentation,  settling,  and  a  slow 
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response  to  the  electric  field.  The  carrier  fluid  is  a  low  viscosity,  low  dielectric 
constant,  high  dielectric  breakdown  strength,  low  conductivity,  nonpolar  fluid. 
Examples  include  kerosene,  mineral  oil,  silicone  (polydimethyl  siloxane), 
greases,  aldehydes,  or  halogenated  hydrocarbons.  Weiss  et  al.  [7]  and  Block  and 
Kelly  [10]  have  provided  a  whole  spectrum  of  different  ER  fluid  formulations. 

The  base  fluid  and  particulate  mixture  also  includes  additives,  namely  surface 
agents,  and  sometimes  an  activating  agent  (water  or  other  polar  solvents  such  as 
alcohols  or  glycols).  One  of  the  roles  of  surfactants  is  to  provide  stability  (to 
deagglomerate  and  to  prevent  sedimentation)  to  the  colloidal  suspension  and 
particle  fluid  interface.  Another  role  is  to  provide  the  particle  surfaces  with 
mobile  charge  carriers  or  an  artificially  large  dielectric  constant. 

The  ER  effect,  either  under  static  or  steady  shearing  conditions,  involves 
organizing  the  suspended  particles  into  multistranded  chains  that  possess 
long-range  order.  The  phenomenon  of  particle  fibrillation,  or  formation  into 
fibrils,  and  subsequent  theory  stemming  from  it,  cannot  fully  account  for  the 
effect  alone  [11].  Specifically,  the  timescale  needed  for  fibrillation  is  significantly 
greater  than  the  typical  millisecond  response  times  of  fluids.  More  recent 
theories  agree  that  the  chaining  is  the  result  of  interparticle  electrostatic 
attraction  that  arises  from  the  particles  polarizing.  Whatever  the  cause  of  the 
polarization,  in  most  fluids  there  is  a  relatively  large  permittivity  difference 
(dielectric  constant)  between  the  suspended  particles  and  the  base  fluid. 
Theories  explaining  the  origin  of  this  polarization  include  an  induced 
double-layer  charge  cloud  surrounding  each  particle  [12,  13],  electro-osmosis- 
activated  water  bridge  [14],  and  intrinsic  particle  polarization  by  electronic, 
atomic,  rotational,  nomadic,  or  interfacial  charge  carriers  [10]. 

The  viscosity  change  is  reversible,  with  an  activation-deactivation  response  time 
of  about  a  few  milliseconds.  When  active,  the  macroscopic  mechanical  and 
electrical  properties  such  as  shear  stress,  viscosity,  and  current  density 
predominantly  depend  on  factors  such  as  electric  field  intensity,  shear  rate, 
temperature,  and  composition  (particle  concentration  and  morphology). 

In  light  of  the  desired  benchmark  properties,  a  multitude  of  component 
parameters  must  be  considered  and  solved  in  the  engineering  and  developing  of 
ER  fluids.  Some  of  the  issues  that  have  seen  some  recent  improvement  include 
(1)  the  compatibility  of  the  base  fluid  and  particulates,  (2)  the  density  and 
dielectric  strength  of  the  base  fluid,  and  (3)  the  presence  of  water  and  other 
impurities  at  the  fluid/particle  interface. 

In  die  past,  the  increasing  demand  for  higher  activity  fluids  was  usually  achieved 
by  higher  particle  concentrations  that  led  to  excessive  agglomeration,  poor 
dispersability,  and  sedimentation.  Associated  with  an  increased  particle 
concentration  is  an  unwanted  increase  of  the  zero-field  viscosity  of  the  fluid. 
Some  of  these  drawbacks  have  been  alleviated  by  better  density  matching 


between  the  components.  Particles  have  become  less  dense,  and  denser-base 
fluids  have  been  introduced.  In  newer  fluids,  silicone  has  been  replaced  by  low 
molecular  weight  (less  than  500)  fluorosilicones  (i.e.,  fluorinated  silicones).  The 
advantages  of  fluorosilicones  are  not  only  based  on  their  higher  densities,  but 
also  on  their  higher  dielectric  strengths  which  permit  higher  electric  fields. 

Since  introducing  anhydrous  fluids  [15,  10],  current  trends  have  taken  a  sharp 
turn  away  from  using  water  to  activate  fluids.  Although  traditionally  water- 
activated  fluids  have  very  high  activities,  using  water  caused  problems  with 
corrosion,  irreproducibility,  limited  temperature  range,  and  instability.  While 
the  amount  of  water  in  "substantially"  anhydrous  fluids  is  still  being  disputed, 
these  fluids  have  made  steady  progress  in  improving  their  activity  levels.  At  the 
same  time,  some  of  these  fluids  tend  to  contain  precursors  (e.g.,  silicates  or 
zeolites)  that  are  highly  abrasive. 

Advancement  has  also  been  made  in  synthesizing,  fabricating,  and  processing 
fluid  components  with  better  control  over  purity,  and  with  particles,  over 
morphology  and  size  distribution.  Introducing  new  types  of  composite  particles 
(e.g.,  metal  coated  polymers)  has  also  resulted  in  improved  density  matching  as 
well  as  reduced  current  densities. 

2.2  Dynamic  Shear  Stress  and  Viscosity  of  the  ER  Fluid 

The  mechanical  behavior  of  the  ER  fluid  may  be  divided  into  a  preyield  and  a 
postyield  regime  [9].  With  increasing  field  intensity,  under  preyield  conditions, 
the  fluid  first  behaves  as  a  viscous,  then  as  a  linear  viscoelastic,  and  finally  as  a 
nonlinear  viscoelastic  solid.  Measurements  of  the  complex  shear  modulus  (i.e., 
storage  and  loss  components)  with  electric  field,  strain  amplitude,  and  frequency 
indicates  that  the  fluid  is  a  relatively  lossy  system.  However,  with  increasing 
electric  field  and  frequency,  the  storage  modulus  correspondingly  increases.  Of 
course,  when  compared  to  more  common  viscoelastic  solids,  the  absolute 
magnitudes  of  the  storage  moduli  are  considerably  less. 

Under  the  postyield  regime,  ER  fluids  are  usually  modeled  as  a  Bingham 
material.  In  steady-state  shear  for  a  Bingham  plastic,  the  shear  stress  is  defined 
as 


x(E,  y)  =  Ty(E)  +  PpY,  (1) 

where  Ty  is  an  electric-field-dependent  yield  stress,  pp  is  plastic  viscosity,  and  y 
is  the  shear  strain  rate.  Typical  experiments  measure  the  shear  stress  as  a 
function  of  shear  strain  rate.  The  intercept  and  slope  correspond  to  ty  and  pp. 
Without  an  electric  field,  there  is  no  ER  effect  and  fluid  behavior  is  Newtonian. 
Some  fluids  may  show  a  small  zero-field  yield  stress,  though.  Thus,  to  a  good 
approximation,  xy  is  zero,  and  pp  is  equal  to  the  Newtonian  viscosity,  p.  This  is 
also  the  viscosity  at  zero-field,  po,  which  is  used  to  characterize  the  fluid. 
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Unlike  the  static  yield  stress  which  varies  linearly  with  the  applied  field 
intensity,  the  dynamic  yield  stress  of  an  ER  fluid  can  be  modeled  by  a  power  law, 
but  usually  a  simple  parabolic  relation  is  sufficient  [5, 16, 17], 

xy  =  AE2.  (2) 

The  proportionality  constant.  A,  quoted  in  [kPa(kVmm-1)'2],  can  be  used  as  the 
fluid  performance  index.  A  value  of  0.25  for  A  implies  a  yield  stress  of  4  kPa  at 
an  electric  field  of  4kVmm‘1.  A  has  been  found  to  vary  with  fluid  and 
temperature. 

While  the  yield  stress  defines  the  upper  limit  of  shear  stress  that  an  ER  fluid  can 
essentially  support,  the  zero-field  viscosity  of  an  ER  fluid  is  important  because  it 
identifies  the  minimum  shear  stress  necessary  to  initiate  flow.  The  temperature 
dependence  of  the  zero-field  viscosity  of  ER  fluids  is  similar  to  that  of  pure  oils 
(i.e.,  it  monotonically  decreases  with  increasing  temperature)  [18].  The 
temperature  dependence  of  the  plastic  viscosity  is  similar,  and  it  varies  with  the 
square  of  the  electric  field.  However,  the  correlations  between  such  relations  are 
not  as  good  as  those  for  the  yield  stress. 

2.3  DC  Current  Density  of  the  ER  Fluid 

Knowing  the  current  density  of  a  fluid  is  critical  because  it  reflects  the  expected 
heat  load  into  the  fluid  and  thereby  affects  its  temperature  stability.  Ultimately, 
the  current  density  determines  the  high  voltage  power  demand  by  the  ER-based 
device.  Because  of  the  inherent  nature  of  the  various  forms  of  polarization 
mechanisms,  conduction  in  ER  fluids  is  non-Ohmic. 

The  measurements  of  current  density  are  fitted  to  quadratic  empirical  relation 

j  =  0.1  (P*E  +  Q*E2),  (3) 

where  j  is  current  density,  E  is  the  field  intensity,  and  P  and  Q  are  constants 
found  by  regression  analysis.  PE  appears  to  be  reminiscent  of  an  Ohmic 
conduction  term.  Origins  of  the  QE2  term  may  be  rooted  in  Coulombic 
interactions  between  particles.  The  role  of  the  0.1  is  rather  arbitrary,  primarily 
for  convenience.  Thus,  a  Q  value  of  10  results  in  16  pAcm  2  at  4  kVmxn  1. 

Q  is  also  function  of  the  temperature  through  an  Arrhenius-type  relation, 

Q  =  Q0exp[-AH(RT)-i],  (4) 

where  R  is  the  Universal  Gas  Constant  and  8.314  J(mol  K)"1  and  AH  represents  an 
activation  energy  over  a  barrier  that  must  be  overcome  in  ER  fluid  conduction. 
For  fluids  that  obey  this  relationship,  logQ  is  a  linear  function  of  the  inverse 
absolute  temperature.  For  water  activated  fluids,  logQ  is  a  linear  function  of  the 
water  content  of  the  particulate  solid  as  well. 

It  is  intuitively  obvious  that  in  a  practical  situation,  an  uncontrolled  ER  fluid 
with  a  large  current  density  may  be  subject  to  a  thermal  runaway  effect.  That  is. 


6 


once  a  fluid  can  no  longer  be  in  thermal  equilibrium  (i.e.,  the  heat  generated 
during  continuous  shearing  cannot  be  transported  out),  the  ensuing  temperature 
rise  will  necessitate  increased  power  levels  to  maintain  the  same  shear  stress 
levels.  In  turn,  this  will  generate  larger  currents. 

Varying  A  with  temperature  can  be  brought  about  indirectly  by  varying  logQ 
with  temperature  (i.e.,  varying  Q  by  changing  the  temperature,  or  varying  Q  by 
changing  the  water  content).  Both  have  the  same  effect  on  A  within  the  limits  of 
experimental  error. 

2.4  Response  Time  of  the  ER  Fluid 

The  activation  response  time  is  linked  to  the  formation,  breaking,  and  reforming 
of  particle  chains.  Conversely,  the  deactivation  response  time  reflects  the  rate  of 
relaxation  of  these  chains  in  the  fluid.  Since  the  strength  of  particle  interaction  is 
determined  by  the  extent  of  polarization,  the  response  time  primarily  depends 
on  the  applied  field.  It  also  depends  on  the  concentration  and  size  of  the 
particles,  as  well  as  the  zero-field  viscosity  of  the  base  fluid  that  affect  the  extent 
of  crowding  and  ease  of  particle  movement. 


3.  Description  of  the  Instrumentation 


The  mechanical  property  measurements  of  ER  fluids  fall  along  one  of  two 
methods:  (1)  sliding  plate  mode  and  (2)  flow  (valve)  mode.  In  the  sliding  plate 
configuration,  the  test  sample  is  constrained  between  a  stationary  and  a  moving 
plate  [16].  Usually  in  the  test,  torque  (and  based  on  the  type  of  fluid  model,  the 
shear  stress)  induced  by  the  sample's  shear  resistance  on  the  stationary  plate  is 
measured  as  a  function  of  the  (constant  or  oscillatory)  speed  of  a  moving  plate 
and  electric  field.  This  type  of  geometry  includes  flat  parallel  plate  arrangements 
and  coaxial  cylinder  tests  either  in  rotary  or  axial  modes.  In  the  flow  mode,  the 
plates  remain  stationary  while  the  test  sample  is  squeezed  through.  Both 
methods  lend  themselves  to  measuring  the  preyield  and  postyield  regions  of  the 
ER  response. 

Response  time  measurement  methods  include  the  rheo-optical  approach, 
oscillatory  parallel  plate  rheometry,  optical-lever/ torque  tube  technique,  and 
pressure-driven  flow  device  [19, 7].  In  all  systems,  either  in  moving  plate  or  flow 
geometry,  the  rate  of  rise  of  shear  stress,  viscosity,  or  pressure  of  the  fluid  is 
matched  to  that  of  the  input  voltage  pulse,  and  the  corresponding  time  delay 
between  signals  is  measured. 
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3.1  ERFD  Constant-Stress  Rheometer 


Conventional  rheometers  use  torsion  springs  to  determine  the  drag  torque 
produced  by  a  fluid.  A  stiff  spring,  which  would  be  adequate  for  measuring  the 
torque  exerted  by  an  "active"  ER  fluid,  would  be  too  insensitive  to  measure  the 
same  fluid  in  its  "inactive"  zero  field  state.  The  relatively  low  shear  rates,  which 
are  needed  for  realistic  measurements  and  achieved  in  such  instruments,  also 
lead  to  excessive  heating  and  cause  the  fluid  to  rapidly  overheat. 

ERFD  has  built  a  Couette  geometry  device  that  does  not  use  torsion  springs  to 
determine  the  exerted  torque.  Owing  to  its  design,  the  measurement  of  the  yield 
stress  is  absolute,  depending  only  on  the  apparatus  dimensions  and  the  applied 
stress.  Similarly,  there  is  no  stick-slip  (problems  at  low  rotational  speeds,  i.e., 
rheopexy),  so  resonances  common  with  torsion  springs  are  also  avoided. 
Viscosities  as  low  as  30  mPa  and  yield  stresses  up  to  7  kPa  can  be  measured. 
Because  the  data  acquisition  per  data  set  is  completed  in  less  than  1  s,  there  are 
little  or  no  heating  problems  associated  with  continually  shearing  the  ER  fluid. 

The  operational  principle  of  the  dynamic  rig  is  extremely  simple.  A  schematic 
diagram  and  photographs  of  the  drop  viscometer  are  shown  in  Figures  l(a-c).  A 
small  (15-30  ml)  fluid  sample  is  loaded  into  a  gap  between  a  cup  and  a  rotor. 
The  source  of  the  constant  stress  is  a  hanging  weight  that  is  attached  with  a 
string  and  wound  around  the  shaft  of  the  rotor.  A  high  voltage  applied  between 
the  cup  and  rotor  generates  an  electric  field,  thereby  modifying  the  fluid 
viscosity.  Thus,  when  the  weight/  string  is  released,  the  rotor  begins  to  spin.  The 
rotation  speed  depends  on  the  net  acceleration  of  the  falling  weight,  which  in 
turn  is  determined  by  the  gravitational  constant  and  the  viscous  drag  in  the  fluid. 
As  the  field  intensity  increases,  the  fluid  viscosity  increases  the  drag  due  to  the 
ER  effect, ,  and  the  rotational  speed  of  the  shaft  is  reduced.  An  optical  wheel  and 
a  photodetector  monitor  the  fall  rate.  The  photodetector  signal  is  fed  through  a 
data  acquisition  system  into  a  computer.  From  the  resultant  displacement  vs. 
time  curves,  a  fitting  routine  determines  the  angular  velocity  and  acceleration, 
and  in  turn  the  shear  stress  and  plastic  viscosity  of  the  fluid.  The  working 
principles  and  calibration  methodology  of  the  dynamic  rig  are  described  in 
Appendix  A. 

Since  the  ER  fluid  is  in  series  with  the  high-voltage  power  supply,  the  DC  current 
can  be  simultaneously  measured.  The  current  drawn  by  the  fluid  is  monitored 
through  the  internal  current  meter  of  the  power  supply.  The  current  density  is 
then  the  ratio  of  the  current  to  the  surface  area  of  the  rotor.  The  temperature 
range  of  the  fluid  can  be  varied  from  room  temperature  to  70  °C  by  using  a 
proportional  heater.  Temperatures  below  ambient  levels  can  be  reached  using  a 
double-walled  cup,  through  which  a  cooling  fluid  (e.g.,  cold  water  or 
alcohol-carbon  dioxide  mixture)  can  be  circulated.  At  the  same  time,  the  shaft 
and  upper  bearings  are  thermally  shielded  and  insulated  from  the  cup  and  rotor 
to  maintain  a  constant  temperature. 


(b)  (C) 

Figure  1.  (a)  Schematic  diagram  of  the  dynamic  drop-viscometer,  (b)  a  photograph  of  the 
device,  and  (c)  a  close-up. 


3.2  Response  Time  Measurements 

The  response  time  tests  were  performed  on  a  counter-rotating,  double-cup  ER, 
clutch-winding  rig  [20,  21].  The  ER  clutch  is  one  among  four  devices  at  the 
University  available  for  response  time  determinations  [20,  21].  The  winding  rig 
is  designed  to  give  a  quick  general  overview  of  the  fluid  characteristics.  The 
loading  and  unloading  time  response  can  be  observed  by  comparing  the  speed  of 
the  winding  system  with  the  speed  of  the  applied  excitation.  The  rate  of 
acceleration  of  the  belt  provides  some  indication  of  the  yield  stress  magnitude. 
The  range  of  shear  rates  is  0  to  about  15,000  s  F  Because  of  the  presence  of  large 
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centrifugal  forces,  the  rig  can  be  used  to  identify  problems  associated  with  fluid 
settling.  A  range  of  temperatures  can  be  covered;  however,  for  temperatures 
above  45  °C,  an  external  blower  is  necessary.  See  Figure  2  for  a  schematic 
diagram  of  the  winding  rig. 


(a) 


Cross  Section  Drawing  of  E,R.  Clutch 


Figure  2.  Schematic  diagram  of  the  double  ER  clutch  (A,  B)  used  to  measure  the  response 
time  is  shown  in  (a),  with  assembly  detail  in  (b).  In  Figure  2(a),  the  various 
components  are  C  -  input  rotor,  D  -  ER  fluid,  E  -  output  rotor,  F  -  pulley,  and 
G  -  drive  belt. 
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The  mechanism  consists  of  two  vertically  mounted  ER  clutches,  A  and  B, 
respectively.  Each  clutch  consists  of  an  input  (C)  and  an  output  (E)  rotor.  The 
gap  between  the  two  rotors  is  filled  with  fluid  (D).  A  pulley  (F),  attached  to  the 
shaft  of  output  rotor  of  clutch  A,  is  connected  by  a  reciprocating  belt  (G)  to  an 
equivalent  pulley  on  the  shaft  of  the  output  rotor  of  clutch  B.  The  motor-driven 
input  rotor  runs  at  virtually  a  constant  speed,  which  is  accomplished  because  its 
inertia  is  much  larger  than  that  of  the  output  rotors,  pulleys,  and  belt.  Because 
the  inertia  of  the  output  rotors,  pulleys,  and  belt  are  as  low  as  possible,  the 
reciprocating  action  of  the  device  can  respond  quickly  to  changes  in  the  ER 
torque  applied  to  the  clutches.  A  full  cycle  is  produced  by  alternatively 
energizing  the  A  or  B  clutch.  Table  2  summarizes  the  various  segments  during 
the  traversing  cycle  of  the  winding  rig. 


Table  2.  Conditions  during  the  traversing  cycle  of  the  winding  rig. 


Portion  of  Belt 
Motion 

Belt  Speed 

Clutch  A 

“  Clutch  B  j 

Energized 

ER  Torque 

Slipping 

Energized 

ER  Torque 

Slipping 

Otol:  constant 
speed 

-u 

no 

no 

yes 

yes 

yes 

no 

1:  excitation  A  on; 
Boff 

-u 

yes 

no 

no 

no 

no 

1  to  2:  slowing 

-u  <  v  <  0 

yes 

no 

yes 

no 

no 

yes 

2:  ER  torque  A  on 

-u  <  v  <  0 

yes 

yes 

no 

no 

yes 

2  to  3:  slowing 

-u<v<0 

yes 

yes 

no 

no 

yes 

3:  stopped 

0 

yes 

yes 

KSH 

no 

no 

yes 

3  to  4:  accelerating 

0  <  v  <  u 

yes 

yes 

yes 

no 

no 

yes 

4  to  6:  constant 
speed 

u 

yes 

no 

no 

no 

yes 

6:  excitation  A  off; 

B  on 

u 

no 

no 

no 

yes 

no 

yes 

6  to  7:  slowing 

no 

no 

yes 

yes 

no 

yes 

7:  ER  torque  B  on 

0  <  v  <  u 

no 

no 

yes 

yes 

yes 

yes 

7  to  8:  slowing 

0  <  v  <  u 

no 

no 

yes 

yes 

yes 

yes 

8:  stopped 

0 

no 

no 

MSM 

yes 

_ _ 

yes 

8  to  9:  accelerating 

-u  <v<0 

no 

no 

yes 

yes 

yes 

yes 

9  to  0:  constant 
speed 

-u 

no 

no 

yes 

yes 

yes 

no 

The  input  rotor  cups  are  motor  driven  to  rotate  in  opposite  directions.  If  the 
viscous  torques  (belt  tension  and  inertial  effects)  are  balanced  on  both  clutches, 
there  is  no  net  torque  on  the  belt,  and  the  output  rotors  are  stationary. 
Otherwise,  the  rotors  spin  either  clockwise  or  counterclockwise.  When  a  step 
voltage  is  applied  to  one  of  the  clutches  (i.e.,  that  ER  clutch  is  activated),  both 
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rotors  rotate  in  that  direction.  If  the  voltage  is  turned  off  and  an  equal  voltage  is 
applied  to  the  other  (inactive)  clutch,  the  rotors  and  belt  rotate  in  the  reverse 
direction.  This  reversal  does  not  occur  instantaneously  because  of  delays  in  the 
response  times  of  the  power  supply,  belt  hydraulics,  and  the  ER  effect.  By 
monitoring  the  speed  of  the  belt  simultaneously  with  the  magnitude  of  the 
applied  voltage  to  either  clutch  as  a  function  of  time,  the  time  delay  of  the  belt 
response  can  be  measured.  Because  the  voltage  rise  time  in  the  current  power 
supply  is  rapid  (i.e.,  less  than  a  millisecond),  most  of  the  delay  could  be 
attributed  to  the  ER  response  (in  the  millisecond  range)  and  thus  be  easily 
determined.  Repeated  tachometer  measurements  and  data  averaging  minimize 
inertial  effects  and  belt  slipping.  The  ER  torque  and  the  electronic-hydraulic 
time  delay  are  both  affected  by  temperature,  shear  rate,  and  field  intensity. 


4.  Experimental  Procedure 


4.1  ERFD  Constant-Stress  Rheometer 

Operating  the  viscometer  was  relatively  straightforward.  Initially,  the  heaters, 
power  supply,  and  computer  were  turned  on.  The  upper  heater,  controlling  the 
temperature  of  the  capstan  shaft,  was  set  at  a  constant  value  of  30  °C,  while  the 
lower  heater,  controlling  the  temperature  of  the  rotor,  holding  cup,  and  fluid, 
could  be  varied  between  30  and  70  °C.  The  fluid  sample  was  pippetted  into  the 
gap  between  the  rotor  and  cup.  After  ensuring  that  no  air  bubbles  remained  in 
die  fluid,  the  temperature  was  set.  Usually,  the  temperature  of  the  cup  and  fluid 
equilibrated  about  30  min  after  its  value  was  set  at  the  proportional  heater.  Six 
temperatures  were  chosen:  5,  15,  30,  40,  50,  and  60  °C.  For  the  temperatures 
below  ambient  conditions,  a  jacketed  cup  was  used  that  permitted  cold  tap  water 
to  circulate. 

Prior  to  a  drop,  the  string,  attached  to  the  pan  holding  a  10-lb  dropweight,  was 
wound  up  around  the  capstan.  A  holding  clamp  was  used  to  lock  the  weight 
into  position.  The  high  voltage  was  turned  on,  and  the  data  logger  was  engaged 
to  receive  data.  By  releasing  the  holding  clamp,  the  weight  was  then  dropped. 
The  total  DC  current  drawn  through  the  fluid  was  recorded,  and  the  high  voltage 
was  turned  off.  During  each  drop,  600  time-pulse  data  points  were  taken  and 
recorded.  A  set  of  six  drops  was  taken  for  each  electric  field  increment.  Between 
the  successive  drops,  the  shaft  was  rotated  60°  counterclockwise,  which  was 
necessary  for  averaging  any  circumferential  nonuniformities  around  the  shaft 
and  bearings  of  the  device. 

The  power  supply  voltage  was  incremented  from  0.0  to  2.0  kV  in  steps  of 
0.25  kV,  which  corresponded  to  a  field  intensity  of  0.5  kVmm’1  at  the  gap.  The 
field  was  ramped  until  either  an  electric  discharge  occurred  in  the  fluid,  or  the 
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drop  weight  did  not  drop.  If  the  weight  did  not  drop,  the  shear  strength  of  the 
fluid  exceeded  the  gravitational  force.  In  those  instances,  the  ramping  sequence 
was  terminated. 

After  completing  each  dataset,  the  computer  was  instructed  to  extract  the  data 
from  the  logger.  Usually,  while  applying  the  high  voltage,  the  computer  was 
disengaged  from  the  logger  to  protect  it  from  any  stray  high-voltage  arcs.  This 
way,  if  the  fluid  arced,  only  the  logger  was  affected.  For  each  dataset,  the 
computer  displayed  the  individual  displacement  vs.  speed  trajectories. 
Trajectories  that  fell  more  than  ±5%  outside  of  the  average  were  typically 
rejected.  A  straight-line  fit  extracted  the  slope  and  intercept  that  could  then  be 
converted  into  shear  stress  and  plastic  viscosity  with  the  appropriate  calibration 
factors.  See  Appendix  A  for  further  details. 

A  new  fluid  sample  was  used  at  each  temperature  increment.  Additionally,  if 
the  fluid  arced,  it  was  removed  from  the  apparatus.  This  was  necessary  because 
a  discharge  across  the  cup/ rotor  gap  usually  damaged  the  rotor  surface  and 
produced  a  hard  lump  in  the  fluid.  With  each  fluid  replacement,  the  cup  and 
rotor  were  rinsed  clean  with  acetone,  then  wiped  clean  with  a  damp  cloth  to 
remove  any  residue.  After  arcing,  the  rotor  was  inspected,  and  the  damaged 
surface  was  repaired  if  necessary. 

The  DC  current  was  recorded  just  prior  to  releasing  the  weight.  The  readings 
from  the  six  consecutive  drops  were  averaged,  converted  into  current  density, 
and  a  straight  line  was  fitted  to  data  to  obtain  the  value  of  Q  at  each  temperature. 
The  current  drawn  by  the  stationary,  nonshearing  fluid  was  usually  less  than 
that  drawn  during  shearing.  This  transient  current  rise  was  attributed  simply  to 
the  rise  of  the  fluid  temperature  from  shearing,  since  a  warmer  fluid  would  draw 
higher  current  levels.  Hence,  the  recorded  current  actually  represented  a  lower 
limit  of  the  electrical  load  across  the  fluid. 

Current  measurements  were  unaffected  by  the  breakdown  of  the  fluid  by  arcing. 
If  arcing  occurred  that  prevented  the  drop  tests  from  continuing,  readings  could 
still  be  taken  as  a  function  of  the  applied  field.  In  these  situations,  the  rotor  was 
spun  around  several  times  by  hand  to  break  up  the  particle  aggregates,  thereby 
rehomogenizing  the  suspension.  In  these  cases,  typically  only  a  single  current 
reading  was  taken.  These  data  were  also  converted  into  current  densities  and 
included  in  the  regression  fit. 

ERFD's  method  of  obtaining  P  and  Q  from  the  equation  for  j  is  to  perform  a  first- 
order,  linear,  least-squares  regression  to  the  empirical  function 

jE-i  =  0.1  (P  +  Q*E).  (5) 

The  linear  term  represents  the  "customary"  Ohmic  conduction  of  the  fluid.  The 
quadratic  term  represents  the  coulombic  interactions  in  the  ER  fluid.  To  a  large 
extent,  P  can  be  treated  as  an  artifact  of  the  fitting  process  [22].  It  is  the  quadratic 


13 


term,  in  which  the  ER  effect  is  related  to  temperature  through  Q,  that  determines 
the  utility  and  working  range  of  the  fluid. 

To  examine  the  effect,  any  artificially  introduced  bias  of  the  fit  near  j  =  0.0  (a 
second-order  least  squares  regression,  wherein  P  and  Q  are  treated 
independently)  was  applied  directly  to  j.  The  results  of  the  two  methods  to 
obtain  P  and  Q  agreed  with  each  other  within  the  limits  of  statistical  error.  The 
Q  values  from  both  methods  were  then  fitted  to  the  exponential  Arrhenius 
temperature  function. 

4.2  Response  Time  Measurements 

The  amount  of  fluid  available  for  testing  was  limited.  Therefore,  only  a  few 
different  conditions  could  be  tried;  Table  3  identifies  these  conditions.  In  all 
cases,  the  signal  frequency  of  the  input  rotors  was  a  2-Hz  binary  square  wave  of 
either  ±1  or  +3kVmm-1.  The  electric  field  intensity  was  based  on  an  input 
rotor/ output  rotor  gap  of  0.5  mm.  In  addition  to  those  listed  in  Table  3,  a  few 
preliminary  tests  were  conducted.  In  these,  the  electric  field  was  applied,  but  the 
rig  was  not  rotating;  or  alternatively,  the  rig  was  rotating  without  any  field 
present. 


Table  3.  Test  matrix  for  the  response  time  measurements. 


Fluid 

T 

rc) 

E 

(kVmm*1) 

Y 

(xlO3  s'1) 

i 

30 

1 

4 

i 

30 

1 

10 

i 

30 

3 

4 

i 

30 

3 

10 

i 

60 

1 

4 

i 

60 

1 

10 

i 

60 

3 

4 

i 

60 

3 

10 

2 

30 

1 

4 

2 

30 

1 

10 

2 

30 

3 

4 

2 

30 

3 

10 

2 

60 

1 

4 

2 

60 

1 

10 

2 

60 

3 

4 

2 

60 

3 

10 

A  ballast  resistor  of  about  50  kQ  was  placed  in  series  with  the  fluid/  electrode  in 
each  clutch  circuit.  The  resistor  did  not  significantly  alter  any  of  the  results.  But 
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it  did  prevent  breakdown  or  arcing  from  occurring,  thereby  giving  a  degree  of 
security  to  the  test. 

An  infrared  thermometer  was  located  on  the  outer  case  of  the  input  rotors.  The 
fluid  temperature  was  assumed  to  be  this  value.  The  ambient  temperature 
ranged  from  20  to  23  °C.  The  belt  speed,  measured  by  a  tachometer,  was  directly 
fed  into  the  digital  oscilloscope.  Because  of  the  characteristic  noise  of  the 
winding  rig  system,  typically  64  readings  of  the  tachometer  output  were  signal 
averaged.  The  on/  off  excitation  response  of  the  two  clutches  was  determined  to 
be  identical.  As  a  result,  only  one  side  of  the  response  was  measured  at  the 
various  field  strengths,  shear  rates,  and  outer  case  temperatures.  Any  imbalance 
in  excitation,  temperature,  or  shear  rate  would  have  caused  the  rig  to  move 
unsymmetrically.  This  was  quite  expected,  since  the  response  time  depends  on 
all  three  variables. 

The  measurement  data  from  the  SMMART  laboratory  was  provided  in  raw 
format  as  a  series  of  graphical  plots  of  traces  of  both  input  and  output  signals  as 
a  function  of  time  for  one  of  the  ER  clutches.  For  each  set  of  conditions 
enumerated  in  Table  3,  several  graphs  were  provided  with  varying  degrees  of 
detail  of  the  voltage  signals.  The  delay  times  were  determined  graphically  by 
measuring  the  elapsed  time  along  the  horizontal  time  scale  from  the  onset  of  the 
sharp  drop  at  the  trailing  edge  of  the  input  voltage  to  the  point  where  the  belt 
began  to  decelerate  (and  accelerate  in  the  reverse  direction).  The  deceleration 
rate  of  the  belt,  which  is  qualitatively  related  to  the  fluid  yield  stress,  was  also 
determined. 


5.  Results  and  Discussion 


The  limited  number  of  currently  available  ER  fluid  formulations  with  their 
unique  physical  characteristics,  such  as  color  and  odor  of  fluid  and/  or  color  of 
suspended  particles,  make  it  relatively  easy  to  correctly  identify  the  supplier  or 
manufacturer  of  a  particular  fluid.  While  these  features  constitute  a  significant 
part  of  a  characterization  effort,  it  would  offer  a  "dead  give-away"  in  guessing  a 
fluid's  composition.  Therefore,  these  aspects  have  been  purposely  omitted  to 
render  the  comparative  measurements  more  unbiased. 

During  viscometry.  Fluid  1  was  fairly  well  behaved.  That  is,  this  fluid  only 
began  to  break  down  (arc  [discharge]  across  the  gap  between  the  rotor  and  cup) 
at  the  highest  temperatures.  However,  Fluid  2  prematurely  arced  at  all 
temperatures,  and  more  importantly,  at  relatively  low  electric  fields.  (An 
interfacial  breakdown  between  the  carrier  oil  and  particulates  in  the  fluid  was 
suspected.  The  actual  reason,  traced  to  contaminants  in  the  oil,  was  only 
revealed  by  GA  and  the  fluid  manufacturer  at  a  much  later  time.) 
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5.1  ERFD  Results 


The  constant-stress  rheometer  was  used  to  measure  the  DC  current  density 
drawn  by  each  fluid,  yield  stress,  and  plastic  viscosity.  The  results  of  the  current 
measurements  will  be  discussed  first,  followed  by  the  stress  and  viscosity 
measurements. 

DC  Current  Density.  Results  of  the  current  density  measurements  for  Fluid  1  are 
summarized  in  Table  4  and  plotted  in  Figures  3-5  as  a  function  of  electric  field 
and  temperature.  Data  for  Fluid  2  are  tabulated  in  a  similar  fashion  as  those  for 
Fluid  1  in  Table  5  and  in  Figures  6-8.  Listed  in  the  tables  are  the  Q  and  P  values 
(Coulombic  and  Ohmic  conductions,  respectively)  from  the  ERFD  and  ARL  least 
squares  fits  to  each  set  of  current  density  data  as  a  function  of  field  intensity  and 
temperature.  In  each  set  of  three  figures,  the  first  two  plots  depict  the  behavior 
of  the  current  density  (j)  and  current  density  per  unit  field  (JE'1)  with  field  (E), 
while  the  third  plot  depicts  the  behavior  with  temperature  (T). 


Table  4.  Fluid  1  current  density  coefficients. 


Temperature 

P 

Q 

ERFD 

ARL 

ERFD 

ARL 

°c 

K 

KD 

(]iAcm'2)(kVmm"1)'1 

(pAcm'2) 

(kVmm-1)'2 

5 

278.16 

0.003,595 

-0.072 

-1.54 

1.13 

1.41 

15 

288.16 

0.003,470 

-1.42 

-3.03 

2.44 

3.00 

30 

303.16 

0.003,299 

-4.03 

-5.96 

9.93 

10.65 

40 

313.16 

0.003,193 

-5.36 

-6.28 

15.59 

15.96 

50 

323.15 

0.003,094 

-17.59 

-43.63 

53.48 

63.59 

60 

333.16 

0.003,002 

-20.04 

-61.70 

106.94 

134.29 

Figures  3  and  4  demonstrate  the  quadratic  increase  of  the  current  with  electric 
field  (i.e.,  applied  voltage)  for  Fluid  1.  Similar  graphs  in  Figures  6  and  7  show 
the  current  for  Fluid  2.  Note  that  true  linear  behavior  only  begins  at  or  above 
1.0  kVmm'1.  These  results  show  two  significant  differences  between  the  fluids. 
First,  below  40  °C,  the  current  density  of  Fluid  1  is  considerably  lower  than  that 
of  Fluid  2.  Second,  despite  the  very  low  currents  at  lower  temperatures,  the  rate 
of  rise  of  the  current  with  temperature  for  Fluid  1  is  significantly  greater  than 
that  of  Fluid  2  (see  Figures  5  and  8  for  a  comparison).  Nevertheless,  it  is  quite 
obvious  that  Fluid  1  would  generate  less  internal  heat  than  Fluid  2  at  or  below 
room  temperature. 
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Fluid  2  Current  Density 


Figure  6.  j  vs.  E;  fluid  2. 


Fluid  2  Current  Density 


9 


Fluid  2  Current  Density 


T  [0CJ 

Figure  8.  j  vs.  T;  fluid  2. 


However,  under  nonequilibrium  conditions  such  as  a  thermal  runaway  scenario. 
Fluid  2  would  heat  up  more  slowly  than  Fluid  1.  This  is  clearly  supported  by 
data  in  Figure  9,  where  the  logQ  is  plotted  against  the  reciprocal  of  the  absolute 
temperature.  The  straight  line  generated  by  a  least  squares  fit  is  included  in  the 
graph.  As  seen  in  the  graph,  both  fluids  fit  well  to  an  Arhenius-type  activation 
process.  The  current  doubling  rates  and  activation  energies  obtained  from  the  P 
and  Q  coefficients  are  listed  in  Table  6.  Note  that  the  activation  energies  of  the 
fluids  lie  between  the  dissociation  energy  of  a  hydrogen  bond  (about 
30  kj/mol'1)  and  that  of  a  "true"  chemical  bond  (approximately  400  kj/ mol'1). 

The  current  doubling  rate  of  8.5  for  Fluid  1  is  slightly  higher  than  the  rate  of 
typical  fluids.  Since  the  doubling  interval  is  unique  to  every  fluid  and  remains 
fairly  constant  with  temperature,  it  appears  that  Fluid  1  was  intended  to  operate 
below  50  °C.  However,  probably  the  most  remarkable  aspect  of  Fluid  2  emerges 
from  Table  6.  Its  current  doubling  rate  is  a  little  more  than  30,  or  roughly  3.5 
times  the  rate  of  Fluid  1.  Thus,  at  higher  temperatures,  resistive  heating  would 
be  significantly  lower  with  its  use.  This  would  greatly  benefit  the  design  of 
ATLAS,  as  lower  parasitic  heating  would  reduce  the  cooling  requirements  of  the 
ER  valve. 

Fluid  Breakdown  Voltage.  Table  7  lists  the  maximum  fields  supported  by  each 
fluid  under  ambient  conditions  (standard  atmospheric  pressure  and  humidity 
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logQ  vs.  T1 


r1  [k1] 


Figure  9.  LogQ  vs.  T'1  for  both  fluids. 


Table  6.  Fluid  current  doubling  rate  and  activation  energy. 


Fluid 

Current  Doubling  Rate  at 

4  kVmm"1  and  300  K 

AH 

kj 

i 

8.5 

63  ±5 

2 

30.8 

15  ±2 

level)  during  testing.  It  may  be  assumed  with  certainty  that  in  cases  where  the 
weight  did  not  drop,  the  maximum  field  supported  by  the  fluid  would  be  at  least 
0.5  kV  or  higher.  With  the  exception  of  the  arcing  at  the  highest  temperature 
(60  °C),  Fluid  1  behaved  well  over  the  entire  temperature  range.  In  contrast. 
Fluid  2  could  not  withstand  the  applied  field  and  broke  down  at  much  lower 
voltages.  This  shortcoming  seriously  compromised  the  course  and  relevance  of 
the  testing. 

The  most  likely  cause  of  arcing  is  a  weak,  low-resistance  path  at  the 
particle/base-fluid  interface  that  allows  an  electric  discharge  to  form  between 
electrodes.  During  testing,  several  explanations  arose  for  the  failure  of  Fluid  2. 
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Table  7.  Fluid  breakdown  voltage. 


mm 

V  breakdown 

HQH 

(kVmm*1) 

Fluid  1 

Fluid  2 

5 

Weight  held  at  3,0 

Arced  at  3.0 

15 

Weight  held  at  3.0 

Arced  at  2.5 

30 

Weight  held  at  4.0 

Arced  at  2.0 

40 

Weight  held  at  4.0 

Arced  at  2.5  and  3.0 

40* 

Not  Applicable 

Arced  at  2.5 

50 

Weight  held  at  3.5 

Arced  at  1.5 

60 

Arced  at  2.0 

Arced  at  1.0 

‘  Degassed  at  room  temperature. 


Approximately  equal  amounts  of  both  fluids  were  air  shipped  by  GA  to  ARL. 
On  arrival,  the  bottle  containing  Fluid  2  was  collapsed  (i.e.,  was  under  partial 
vacuum).  Presumably,  a  significant  fraction  of  the  air  was  inadvertently 
adsorbed  onto  the  particle  surfaces  during  shipment  from  San  Diego,  CA,  to 
Aberdeen  Proving  Ground,  MD.  The  fluid  was  degassed  for  a  few  hours  with  a 
single-stage  rotary  pump  to  extract  most  of  the  dissolved  air.  Since  this  pump 
would  likely  not  remove  dissociated  gases,  this  method  was  treated  with 
skepticism.  As  such,  during  the  evacuation  with  the  roughing  pump,  only  a  few 
air  bubbles  emerged  from  the  fluid.  Supported  by  the  lack  of  improvement  in 
the  current  data,  the  effect  of  degassing  resulted  in  a  very  slight  change  of  the 
rheological  properties. 

An  equally  harmful  pathogen  with  deleterious  effects  for  the  fluid  would  have 
been  moisture  from  air.  For  this,  a  vacuum  bakeout  would  have  sufficed; 
however,  ERFD  lacked  the  facility  to  perform  this  operation.  Alternatively,  it 
was  put  forth  that  perhaps,  triggered  by  various  changes  in  ambient  conditions 
(temperature,  pressure,  and  gaseous  contaminants),  a  tertiary  component  in  the 
fluid  (such  as  an  additive)  or  surfactant  reacted  with  the  base  oil,  forming  a  more 
conductive  product  layer  on  the  particles.  Lastly,  it  was  hypothesized  that 
fundamental  problems  already  existed  with  the  fluid,  but  because  of  inadequate 
or  incomplete  testing  prior  to  shipment,  the  breakdown  went  unnoticed. 

Yield  Stress  and  Viscosity.  The  yield  stress  data  as  a  function  of  electric  field  for 
Fluid  1  and  Fluid  2  are  plotted  in  Figures  10  and  11;  the  data  are  plotted  against 
the  square  of  the  electric  field.  As  is  apparent  in  die  figures,  the  yield  stress  of 
Fluid  1  adheres  to  the  quadratic  dependence  fairly  well,  regardless  of  the 
temperature.  In  contrast,  adherence  to  the  expected  behavior  of  Fluid  2  is 
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considerably  poorer,  most  likely  attributed  to  its  premature  breakdown.  The 
deviation  from  the  expected  straight-line  fit  in  the  latter  data  is  greater,  especially 
at  higher  temperatures.  The  arcing  and  subsequent  nonuniform  shearing  of  the 
arced  fluid  causes  bad  drop  trajectories  and  poor  fits,  which  could  easily 
correspond  to  unrealistic  negative  yield  stress  values. 

The  results  in  Figure  11  also  reveal  that  the  rheological  properties  of  the 
outgassed  fluid  did  not  significantly  change  from  those  of  the  as-received  fluid. 
Similarly,  based  on  the  breakdown  voltage  shown  in  Table  7,  the  outgassing 
treatment  did  not  prevent  or  delay  arcing  from  recurring  either. 

Whenever  arcing  occurred,  the  damaged  fluid  was  replaced  with  a  new  batch  of 
fresh  fluid.  Consequently,  there  is  considerably  more  data  for  Fluid  2  at  lower 
electric  fields.  Despite  the  arcing  problems,  the  multiple  data  values  for  Fluid  2 
reflect  fairly  good  measurement  reproducibility. 

With  increasing  temperature  (up  to  60  °C),  the  yield  stress  of  Fluid  1  increased 
slightly.  This  trend  is  clearly  seen  in  Figure  12.  As  revealed  by  the  slope  of  the 
lines,  at  lower  field  intensities  the  increase  is  more  gradual.  At  higher  fields,  the 
increase  is  more  pronounced.  A  similar  trend  is  shown  in  Figure  13  for  Fluid  2. 
But,  for  Fluid  2,  only  data  from  0  to  30  °C  should  be  considered  in  the 
comparison.  Because  arcing  occurs  above  this  temperature,  the  data  is 
incomplete  and  unreliable;  therefore,  only  limited  conclusions  on  the  yield  stress 
should  be  drawn.  However,  over  this  restricted  range,  the  rate  of  increase  of  the 
yield  stress  for  Fluid  2  appears  to  be  greater  than  that  of  Fluid  1. 

The  proportionality  constant,  or  the  strength  factor.  A,  relating  the  electric  field 
to  the  yield  stress,  was  obtained  from  linear  least  squares  fits  to  the  data  plotted 
in  Figures  10  and  11.  Generally,  as  shown  in  Tables  8  and  9  and  Figure  14,  the 
value  of  A  increases  with  temperature  for  both  fluids.  In  Table  9,  the  negative  A 
values  are  only  included  for  the  sake  of  completeness  (i.e.,  those  values  that  arose 
from  fitting  to  negative  stress  values),  but  they  were  excluded  from  the  figure. 
From  0  to  60  °C,  A  for  Fluid  1  increases  roughly  twofold  from  0.15  to  about  0.29. 
It  is  expected  that  A  will  not  increase  indefinitely,  but  will  reach  some 
asymptotic  value.  A  similar  trend  was  observed  for  Fluid  2,  but  its  strength 
factor  was  considerably  higher,  reaching  values  above  0.35.  However,  for  this 
fluid  the  electric  field  could  not  be  increased  above  2.5  kVmm'1  (see  Table  7).  As 
a  result,  the  Fluid  2  fits  used  to  obtain  A  are  based  on  fewer  points;  typically, 
standard  errors  were  about  10-15%,  or  3-5  times  larger  than  those  for  Fluid  1. 

The  yield  stress  at  an  electric  field  of  4kVmm"1  is  a  critical  benchmark  for 
evaluating  any  ER  fluid  candidate  for  recoil  applications.  However,  no  actual 
yield  stress  data  at  4kVmm’1  was  taken  for  either  fluid.  During  the  attempted 
drops.  Fluid  1  was  well  behaved  and  did  not  arc.  In  this  case,  no  drop  data  could 
be  taken  because  the  dropweight  did  not  fall;  that  is,  the  ER  force  was  greater 
than  the  force  of  gravity.  It  is  quite  reasonable  to  assume  that  working  with 


24 


Tmax  [kPa]  Tmax  [kPa] 


Fluid  1  Yield  Stress 


Table  8.  Fluid  1  strength  factor,  calculated  yield  stress  and  no-field  viscosity. 


T 

(°Q 

Emax 

kVnun’1 

A 

kPafkVmin'1)"2 

Xy(2  kVmm'1) 
kPa 

Ty(3  kVmin'1) 
kPa 

xy(4  kVmrrf1) 
kPa 

no 

mPa*s 

5 

3.0 

0.15  ±  0.01 

0.59 

1.3 

2.4 

350 

15 

3.0 

0.14  ±0.02 

0.55 

1.2 

2.2 

270 

30 

4.0 

0.21  ±0.01 

0.85 

1.9 

3.4 

130 

40 

4.0 

0.20  ±  0.01 

0.80 

1.8 

3.2 

-6 

50 

3.5 

0.24  ±0.01 

0.98 

2.2 

78 

60 

2.0 

0.29  ±0.04 

1.16 

2.6 

4.6 

9 

Note:  shaded  areas  correspond  to  extrapolated  values  beyond  measurement  range. 


Table  9 .  Fluid  2  strength  factor,  calculated  yield  stress  and  no-field  viscosity. 


Emax 

kVrrtm"1 

A 

kPatkVmm'1)'2 

xy(2  kVmm"1) 
kPa 

xy(3  kVmm‘>) 
kPa 

Ty(4  kVmm"1) 
kPa 

no 

mPa«s 

5 

3.0 

0.22  ±0.01 

0.89 

2.0 

3.6 

60 

15 

2.5 

0.31  ±0.03 

1.26 

2.8 

5.0 

78 

30 

.0 

0.35  ±  0.03 

1.41 

3.2 

5.7 

85 

40 

.5 

0.20  ±  0.03 

0.80 

1.8 

3.2 

38 

40a 

2.5 

0.30  ±  0.04 

1.20 

2.7 

4.8 

31 

50 

.5 

-0.96  ±1.20 

- 

- 

- 

22 

60 

1.0 

-0.51 

— 

— 

-35 

Note:  shaded  areas  correspond  to  extrapolated  values  beyond  the  actual  measurement  range. 
a  Degassed  at  room  temperature. 


larger  dropweights  would  have  overcome  this  minor  problem.  Unfortunately, 
due  to  time  considerations,  the  additional  tests  with  larger  dropweights  were  not 
performed.  Instead,  the  yield  strengths  were  calculated  using  equation  2.  The 
calculated  yield  strength  of  3.4  kPa  at  4  kVmrrf1  and  30  °C  falls  somewhat  short 
of  the  4-kPa  level. 

No  such  conclusions  could  be  drawn  for  Fluid  2.  The  unavailability  of  high  field 
data,  coupled  with  premature  breakdown,  seriously  compromised  the  extent  of 
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Figure  14.  A  vs.  T'1  for  both  fluids. 


the  comparative  testing.  Consequently,  the  validity  of  the  calculated  yield  stress 
for  Fluid  2  above  2  kVmm"1  is  highly  controversial.  Nevertheless,  based  on  the 
available  low  field  data,  the  yield  stresses  at  3  and  4  kVmm"1  were  calculated  to 
extrapolate  fluid  behavior.  With  that  in  mind,  the  projected  yield  strength  of 
Fluid  2  at  4  kVmm"1  and  30  °C  was  estimated  at  about  4.0  ±1.0  kPa,  which  would 
be  within  the  desired  target  range  (see  Table  9). 

For  a  Bingham  fluid,  the  plastic  viscosity,  qp,  represents  the  slope  of  the  family  of 
shear  strain  rate  curves  (straight  lines)  vs.  the  shear  stress,  and  it  decreases  with 
electric  field  and  temperature.  Figures  15  and  16  illustrate  this  behavior  quite 
well  for  both  fluids.  For  Fluid  1,  the  range  of  measured  viscosities  is  very  wide. 
For  Fluid  2,  the  viscosity  range  is  considerably  narrower.  In  both  cases,  the 
dependence  %  on  the  square  of  the  electric  field  is  generally  linear.  As  depicted 
in  Figure  15,  the  proportionality  constant,  or  the  rate  of  change  of  qp,  decreases  as 
the  temperature  increases.  This  effect  is  not  as  clearly  observable  for  Fluid  2  (see 
Figure  16).  The  anomalous  plastic  viscosity  datum  point  at  50  °C  was  most  likely 
caused  by  arcing  during  the  measurement  which  resulted  in  a  bad  fit  to  the 
trajectory.  It  was  included  in  the  graph  for  the  sake  of  completeness. 

In  practical  terms,  the  zero-field  viscosity,  qo  (i.e.,  qp  without  an  applied  electric 
field),  reflects  the  ease  at  which  the  fluid  can  be  pumped  through  the  hydraulic 
lines  of  the  recoil  device.  At  a  more  fundamental  level,  the  magnitude  of  qo 
depends  in  part  on  the  morphology  and  the  agglomeration  tendency  of  the 
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Fluid  1  Plastic  Viscosity 


Figure  15.  r|p  vs.  E2;  fluid  1. 


Fluid  2  Plastic  Viscosity 


Figure  16.  tip  vs.  E2;  fluid  2. 


particles.  However,  the  base  viscosity  of  the  carrier  oil,  primarily  determined  by 
its  molecular  weight  and  structure,  has  a  more  significant  impact  on  r)o.  As 
expected,  Tables  8  and  9  indicate  that  the  two  fluids  are  very  different.  The 
experimental  errors  on  the  viscosity  measurements  are  about  ±10  mPa*s. 
Determining  r]o  is  not  subject  to  the  issues  associated  with  premature  breakdown 
caused  by  arcing.  The  zero  or  negative  values  of  r)o  are  artifacts  of  the  fitting 
process  and  would  most  likely  be  below  the  measurement  limit  of  the 
instrument.  As  shown  in  Figure  17,  at  0  °C,  the  viscosity  of  Fluid  1  is  about 
350  mPa*s,  while  at  60  °C  it  is  near  9  mPa*s.  In  contrast,  the  viscosity  of  Fluid  2 
is  about  60  mPa»s  at  0  °C  and  less  than  20  mPa»s  at  60  °C.  At  the  lower  end  of 
the  temperature  scale.  Fluid  1  is  much  thicker  than  Fluid  2.  Based  on  its  low 
temperature  viscosity,  it  is  suspected  that  Fluid  1  would  be  rather  difficult  to 
circulate  below  0  °C;  however,  at  higher  temperatures,  circulation  would  not  be  a 
problem.  In  contrast,  the  rjo  for  Fluid  2  is  not  only  lower,  but  also  it  does  not  vary 
as  much  as  that  of  Fluid  1.  In  turn,  the  lower  temperature  dependence  of  the 
fluid  viscosity  allows  the  use  of  a  less  robust  pumping  system. 
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Figure  17.  r)o  vs.  T'1  for  both  fluids. 

Based  on  the  current  draw,  zero-field  viscosity,  and  extrapolated  shear 
properties,  it  may  be  concluded  that  once  corrected.  Fluid  2  will  be  closer  to 
meeting  the  benchmark  conditions  than  Fluid  1.  But  it  appears  that  Fluid  1  was 
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designed  to  operate  within  a  narrow  temperature  regime,  between  30  and  50  °C. 
Below  this  range,  the  zero-field  viscosity  is  too  high,  whereas  above  this  range, 
the  current  draw  is  too  high.  Provided  adequate  cooling  is  available,  this  range 
could  be  extended  to  higher  temperatures.  However,  using  Fluid  1  at  lower 
temperatures  will  be  greatly  limited  by  the  rapid  increase  in  its  zero-field 
viscosity. 

5.2  University  of  Sheffield  Results 

Concurrent  with  the  testing  at  ERFD,  samples  of  both  fluids  were  taken  to  the 
SMMART  laboratory  at  the  University  of  Sheffield.  Out  of  courtesy,  the 
laboratory  was  informed  of  the  premature  arcing  encountered  while  testing 
Fluid  2.  Specifically,  as  a  further  precaution,  the  electric  field  was  to  be  ramped 
up  slowly  to  the  preselected  value  (e.g.,  3  kVmm'1)  to  protect  the  apparatus. 
Upon  receiving  this  recommendation,  researchers  at  the  SMMART  laboratory 
further  consulted  with  Professor  Harry  Block  of  Cranfield  University,  a  known 
authority  on  ER  fluids.  At  Block's  suggestion,  a  vacuum  bakeout  at 
80  °C  was  performed  on  Fluid  2,  which  prevented  arcing  during  subsequent 
winding  rig  tests. 

This  rather  simple  solution  implied  that  the  cause  of  the  breakdown  was  most 
likely  a  low  boiling  point  impurity  or  contaminant  in  the  fluid  (i.e.,  a  low 
molecular  weight  hydrocarbon  or  water).  This  was  also  consistent  with  the  fact 
that  the  impurity  could  not  be  removed  with  a  roughing  pump  at  room 
temperature,  but  was  effectively  driven  off  at  intermediate  temperatures  under 
better  vacuum. 

The  precautionary  step  provided  a  viable  remedy  for  the  eliminating/ reducing 
the  arcing.  However,  there  was  insufficient  fluid  or  time  remaining  to  repeat  the 
entire  range  of  viscometry  experiments.  More  importantly,  this  action  essentially 
nullified  the  initial  purpose  of  the  measurements — comparative  testing  of  the  ER 
fluids  in  their  as-received  states.  Preliminary  experiments  demonstrated  that 
there  was  little  or  no  interference  in  the  measurements.  These  tests  also 
indicated  that  the  measurements  were  repeatable. 

Results  of  the  response  time  measurements  are  listed  in  Table  10.  A  typical 
example  of  the  oscilloscope  traces  is  also  shown  in  Figures  18(a)  and  18(b).  In  the 
figure,  the  upper  trace  is  the  DC  voltage  applied  to  one  of  the  ER  clutches.  The 
lower  trace  is  the  tachometer  output  corresponding  to  the  movement  of  the  belt. 
Recall  that  when  the  applied  field  on  one  of  the  clutches  goes  to  zero,  the  applied 
field  on  the  other  clutch  increases.  As  shown  in  the  lower  trace,  the  belt  is  either 
moving  clockwise  (CW)  (defined  as  positive)  or  counterclockwise  (CCW).  The 
portion  of  the  lower  trace,  where  the  tachometer  output  linearly  increases  or 
decreases,  reflects  the  deceleration  and  acceleration  of  the  belt  associated  with 
reversing  the  direction  of  motion.  The  alternating  belt  rotation  follows  switching 
the  DC  step  voltage  (or  electric  field)  input  with  a  slight  time  delay.  On  the 


30 


Table  10.  Response  time  measurement  results. 


Fluid 

T 

(°C) 

E 

(kVmm'1) 

Y 

(xlO3  s'1) 

tr 

(ms) 

Deceleration  Rate 

(Vs'1) 

i 

32.0 

1 

4 

2.4 

34 

i 

32.5 

1 

10 

3.4 

23 

1 

30.0 

3 

4 

2.6 

140 

i 

34.0 

3 

10 

2.8 

240 

i 

60.0 

1 

4 

2.5 

22 

i 

59.5 

1 

10 

2.5 

34 

i 

62.0 

3 

4 

2.7 

45 

i 

60.0 

3 

10 

2.0 

270 

2 

30.0 

1 

4 

2.0 

35 

2 

32.0 

1 

10 

1.2 

5 

2 

31.0 

3 

4 

2.0 

180 

2 

33.5 

3 

10 

2.1 

140 

2 

60.0 

1 

4 

2.0 

39 

2 

59.0 

1 

10 

2.4 

39 

2 

62.0 

3 

4 

1.6 

200 

2 

60.0 

3 

10 

1.2 

190 

expanded  scale  in  Figure  18(b),  this  time  delay  is  clearly  apparent.  The  DC  step 
voltage  input  usually  drops  from  its  constant  level  to  OkV  within  0.6-1. 0  ms. 
The  delay  time,  tr,  or  ER  response  time,  defined  as  the  elapsed  time  from  the 
onset  of  the  drop  of  the  DC  field  to  the  intersection  between  the  linear  rise  and 
the  horizontal  portion  of  the  belt  speed  trace,  was  not  corrected  for  the  DC 
voltage  decay. 

The  response  time  of  Fluid  1  varies  from  2.0  to  3.4  ms.  In  contrast,  the  response 
time  of  Fluid  2  is  somewhat  less,  varying  from  1.2  to  2.4  ms.  The  response  time 
also  slightly  decreases  for  both  fluids  with  temperature.  However,  with  respect 
to  the  applied  electric  field  or  shear  rate,  there  is  no  clear  trend.  Assuming  that 
all  other  factors  were  the  same,  the  deceleration  rates  listed  in  the  table  reflect 
that  there  is  about  a  four  to  six  fold  increase  in  the  yield  stress  of  the  fluids  from 
1  to  3  kVmm’1.  This  generally  agrees  with  the  constant  stress  viscometer  results, 
which  indicated  that  the  yield  stress  increase  from  1  to  3  kVmm'1  was  roughly  a 
factor  of  ten.  Because  the  winding  rig  only  gave  a  quick,  general  overview  of  the 
fluid  characteristics,  the  data  was  not  analyzed  or  calibrated  further. 
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Figure  18.  Response  time  measured  data  for  T  =  30  °C,  E  =  3kVmm  l,  and 
Y  =  4  x  103  s~l .  DC  voltage  input  and  tachometer  output  are  shown  in 
(a)  and  expanded  in  view  of  tachometer  output  in  (b)  reveals  time  delay 
(attributed  partly  to  the  ER  fluid).  Note:  traces  in  (b)  correspond  to  area 
highlighted  by  the  green  rectangle  in  (a). 
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6.  Recommendations  and  Conclusions 


Results  showed  that  Fluid  1  remained  stable  over  the  testing  range  of  electric 
fields  and  temperatures.  Its  activity  was  near  the  benchmark  criteria.  At  lower 
temperatures,  the  current  density  was  below  or  near  the  desired  target  value.  At 
higher  temperatures,  above  50  °C,  the  current  density  quickly  exceeded  the 
target  value  by  factors  of  five  or  six.  The  data  also  revealed  that  the  fluid  has  an 
acceptable  zero-field  viscosity  range  above  30  °C.  However,  because  of  its  fairly 
high  zero-field  viscosity  below  this  temperature,  future  improvements  will  be 
necessary,  either  by  reducing  the  volume  fraction  of  the  particles  (at  the  expense 
of  activity),  or  by  using  a  lighter,  less  viscous  base  fluid. 

Because  of  the  premature  breakdown  at  relatively  low  electric  fields,  it  was 
difficult  to  realistically  assess  Fluid  2.  Based  on  the  partial  data  at  low  electric 
fields  and/or  extrapolated  results  beyond  the  breakdown  conditions,  the 
following  conclusions  may  be  drawn.  At  low  temperatures,  the  current  density 
was  about  60%  above  the  desired  value.  At  higher  temperatures  and  high  fields, 
the  absolute  value  of  the  current  density  missed  the  target  value  by  a  factor  of 
two.  However,  the  rate  of  increase  of  the  current  density  of  this  fluid  (which  is 
presumed  to  be  independent  of  the  fluid  breakdown),  with  respect  to 
temperature,  was  significantly  lower  than  that  of  Fluid  1.  Fluid  2  met  the 
benchmark  criteria  in  another  category,  the  zero-field  viscosity.  Finally, 
extrapolating  the  high-field  activities  from  those  at  low  fields  indicated  that 
high-field  activities  of  a  "corrected"  fluid  (i.e.,  one  that  does  not  break  down  at 
intermediate  voltages)  would  also  meet  the  benchmark  criteria. 

Overall  response  time  measurement  results  indicated  that  Fluid  2  was  slightly 
faster  than  Fluid  1.  Within  the  limited  range  of  test  variables,  both  fluids  had 
activation/ deactivation  response  times  of  the  order  of  2.5-3.0  ms  for  Fluid  1,  and 
about  2.0  ms  for  Fluid  2.  No  consistent  trend  of  the  response  time  data  was 
observed  for  either  fluid,  with  respect  to  a  temperature,  shear  rate,  or  field 
strength  variation. 

Four  fundamental  yet  undemonstrated  aspects  of  these  fluids  remain,  and  they 
all  have  a  direct  impact  on  implementing  the  ATLAS  program. 

(1)  The  premature  settling  and  sedimentation  caused  by  the  remaining 
density  differences  between  the  particulates  and  the  base  fluid, 

(2)  particle  abrasion  on  hydraulic  seal  integrity, 

(3)  excessive  current  drawn  resulting  in  internal  overheating  and 
subsequent  dielectric  breakdown,  and 

(4)  extendibility  of  the  minimum  service  temperature  to  well  below  0  °C. 
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Both  fluids  must  be  further  tested  to  answer  these  unresolved  issues.  Future 
measurements  of  the  fluid  properties  must  also  consider  the 
changing/unchanging  properties  with  repeated  cycling  and  over  an  extended 
time  period.  Additionally,  the  fluids'  compatibility  with  certain  metals,  alloys, 
and  polymer-based  materials  (e.g.,  copper,  titanium  alloys,  stainless  steel,  brass, 
or  carbon-epoxy  composites)  have  yet  to  be  addressed. 


7.  Postscript 


Since  the  writing  of  this  report,  the  development  of  improved  fluids  has 
continued,  and  newer  formulations/ updates  from  both  manufacturers  have 
become  available.  Based  on  the  manufacturer's  specification  sheets,  fluid 
properties  are  beginning  to  reach  the  desired  goals  set  forth  in  Table  1. 

The  manufacturer  of  Fluid  2  has  released  two  updates  of  the  same  fluid 
formulation.  The  first  update,  released  in  1997,  was  primarily  in  response  to 
solving  and  eliminating  the  cause  of  the  premature  breakdown.  The  cause  of  the 
breakdown  was  related  to  impurities  on  the  particles  and  in  the  base  fluid  when 
a  switch  of  material  suppliers  occurred.  The  properties  of  this  update  were 
identical  to  Fluid  2  in  most  respects,  though  the  activity  level  was  slightly  lower 
than  anticipated;  the  yield  strength  was  about  3  kPa,  or  75%  of  the  desired  value. 
A  new  release  became  available  in  1998,  and  research  has  continued  through 
1999.  The  new  fluid  has  an  acceptable  activity  level  (i.e.,  near  4  kPa),  but  with  a 
somewhat  higher  no-field  viscosity  at  lower  temperatures. 

The  other  fluid  manufacturer  has  also  produced  a  new  experimental  fluid  that 
supposedly  has  the  same  characteristics  as  Fluid  1,  but  with  an  activity  level  that 
is  double  the  original  formulation.  All  three  fluid  samples  await  property 
verification  at  ARL. 
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Appendix  A.  Theory,  Governing  Relations,  and  Calibration 
Method  for  the  Constant-Stress  Rheometer 


The  force  balance  for  the  dropweight  in  the  dynamic  rig  under  the  combined 
action  of  gravity,  drag  of  the  bearings,  and  the  viscous  force  exerted  by  any  fluid 
between  the  cup  and  bob  is  given  as: 

Ma  =  Mg  -  IR‘2a  -  k'v  -  (viscous  torque)R_1. 


Newtonian  Fluid.  For  a  Newtonian  liquid  in  the  cup,  the  viscous  force  is 
defined  as 

F  VISCOUS  =  4nhpv[R2(ri2  -  r22)p  =  Dpv, 

where  p  is  the  viscosity,  R  is  the  capstan  radius,  ri  is  the  radius  of  the  bob,  12  is 
the  radius  of  the  cup,  h  is  the  bob  height,  and  k'v  is  the  drag  term  for  the  bearing. 
Upon  insertion,  the  relation  becomes 

Ma  =  Mg  -  IR‘2a  -  k'v  -  Dpv, 


or 


a  -  Mg[M  +  IR‘2p  -  (k'  +  Dp)[M  +  IR'2p  v. 
The  relationship  of  a  vs.  v  is  linear,  with  an  intercept  of 

ao  =  Mg[M  +  IR~2p 


and  a  slope  of 

k  =  -  (k'  +  Dp)[M  +  IR“2p. 

The  parameters,  called  a  and  k,  are  calculated  by  the  ERFD's  "Lindafit"  data 
reduction  program.  Lindafit  obtains  these  values  directly  from  the  displacement 
vs.  time  trajectory. 

Bingham  Fluid.  For  a  Bingham  plastic,  the  shear  stress  is  defined  as 

T  =  Ty+Jlpy, 

where  xy  is  yield  stress,  pp  is  plastic  viscosity,  and  y  is  the  shear  strain  rate.  In 
the  absence  of  an  electric  field,  xy  is  zero,  and  pp  is  equal  to  the  Newtonian 
viscosity  p.  The  derivation  for  the  dynamic  viscometer  assumes  that  the  fluid  is 
partially  shearing.  That  is,  only  the  portion  of  fluid  is  shearing  where  the  yield 
stress  has  been  exceeded  in  the  gap.  When  this  is  inserted  into  the  basic 
equation,  the  presence  of  a  yield  stress  decreases  the  measured  value  of  a  by  an 
amount  equal  to  Aa.  The  viscosity  term  remains  unchanged  and  acts  through  the 
viscous  torque  term. 
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It  can  be  easily  shown  that  for  the  Bingham  fluid 

Aa  =  4nhn2[2R(M  +  IR‘2)piy. 

Textbooks  show  that  on  non-Newtonian  fluids  (i.e.,  for  conventional  and 
Bingham  plastics),  when  the  angular  velocity  is 

£2  =  Typp'1[r22(2ri)2  -  0.5  -  ln(r2rf  *)], 

all  of  the  material  is  being  sheared.  The  term  in  brackets  is  0.0004  for  ERFD's 
system.  At  xy  =  3  kPa,  pp  =  150  mPa*s,  which  leads  to  a  pulse  length  of  8,000  ps. 
The  first  two  time-interval  pulses  usually  exceed  this  on  the  dynamic  rig,  but 
after  that  the  pulses  are  considerably  less  (i.e.,  all  of  the  material  is  shearing  after 
the  first  1-3  pulses). 

Calibration  Procedure.  The  calibration  procedure  is  performed  using  a  series  of 
drop-weights  and  different  viscosity  oils  with  Newtonian  behavior.  Returning  to 
the  previous  equations  for  Newtonian  fluids,  note  that  the  intercept,  ao,  is  a 
function  only  of  the  weight  used,  not  of  the  temperature  and  not  of  the  cup 
contents.  This  has  been  verified  by  ERFD.  Graphs  of  ao 1  vs.  M"1  are  also  linear, 
and  the  inertial  term  can  easily  be  obtained.  An  absolute  calibration  may  be 
obtained  by  verifying  that  the  inverse  of  the  intercept  of  tius  plot  does 
correspond  to  the  value  of  the  gravitational  constant,  g. 


The  slopes  of  the  calibration  lines  are 

k  =  -  k'[M  +  IR'2p  -  Dp[M  +  IR‘2p. 

The  k  values  depend  on  the  weight  used,  as  well  as  on  the  temperature  and  type 
of  oil  in  the  cup.  (The  dependence  is  through  the  fluid  viscosity,  which  is  a 
function  of  temperature.)  The  remaining  unknown  terms  in  this  equation  can  be 
determined  using  a  series  of  oils  of  known  viscosity.  Based  on  the  previous 
equation,  a  graph  of  k  (as  calculated  by  Lindafit)  against  p  leads  to  a  series  of 
straight  lines  with  an  intercept  of 

k'[M  +  IR‘2p 


and  slope  of 


D[M  +  IR'2p, 


where 


D  =  4nh[R2(ri2-  r22)p. 

Values  for  the  relation  of  k  vs.  p  may  be  refined  using  these  relationships,  but  are 
basically  drawn  directly  from  the  calibration.  In  the  previous  equations, 

h  =  bob  height  =  50  mm, 

ri  =  bob  radius  =  25  mm. 


r2  =  cup  radius  =  25.5  mm. 
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R  =  capstan  radius  =  17.5  mm,  and 

M  +  IR’2  is  already  determined  during  calibration  with  respect  to  the 
dropweight. 

The  best  relation  for  a  10-lb  dropweight  from  the  last  calibration  of  the 
instrument  (performed  January  1995)  is 

p  =  (k  -  0.419)(0.00375)'1, 

and  these  constants  are  temperature  independent.  The  estimated  error  of  p  is 
about  ±10  mPa»s  for  oils  of  known  viscosity. 

As  printed  from  the  Lindafit  calculation,  Aa  is  angular  acceleration,  and  must 
therefore  be  multiplied  by  the  capstan  radius  (0.0175  m)  to  give  units  of  ms‘2. 
Working  out  the  relationship  for  a  Bingham  fluid, 

xy[kPa]  =  10.53  x  10*3  Aa. 

It  is  also  clear  that  when  a  is  0  (or  Aa  equals  the  zero  field  value  of  a),  the  weight 
will  not  fall.  By  calculation,  this  implies  a  yield  stress  of  3.8  kPa.  In  practice, 
because  of  the  difference  between  static  and  dynamic  yield  stress,  runs  with  a 
10-lb  weight  cannot  measure  more  than  3.2-3.3  kPa. 

Using  18-lb  weights  will  allow  measurements  of  up  to  7  kPa.  The  expected  value 
of  a  equals  434  ±  10;  the  calculated  yield  stress  factor  equals  16.2  x  10'3. 
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Table  B-3.  Fluid  2  current  draw 


Run  Temperature 


Field 

(kVmm'i) 


84  83  87 


177  179  177 


282  283  285 


433  424  435 


592  590  602 


0  0  0 


32  32  33 


110  108  110 


222  225  228 


375  372  378 


552  565  552 


0  0  0 


55  57  57 


189  187  186 


368  372  396 


616  615  626 


0  0  0 


71  75  75 


241  250  251 


476  483  486 


790  794  765 


1121  1119 


1460  1510  1440 


1070  1080  1100 


1460  1500  1460 


0  0  0 


74  75  76 


249  254  256 


488  490  495 


785  777  805 


1100  1120  1110 


1530 


25 


87  88 


182  181 


298  298 


454  446 


0  0 


34  32 


113 


232  236 


385  389 


583  575 


0  0 


55  56 


192  196 


377  387 


656 


0  0 


77  75 


246  250 


498  480 


804  777 


0  0 


75  75 


247  I  245 


508 


818  815 


2540 


0  0  0  I  0  I  0  I  0 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Deeassed 


Table  B-3.  Fluid  2  current  draw  (continued) 


Table  B -4.  Fluid  2  yield  stress. 
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